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Abstract 

The functional derivative of the effective action with respect to an external field is 
part of the equation of motion of this field if one-loop effects induced by quantum 
fluctuations or thermal fluctuations are included in minimizing the action of this 
field. Examples occur in all field theories displaying classical solutions or - as 
the Nambu-Jona-Lasinio model - selfconsistent field configurations. We describe 
here a numerical method for computing such functional derivatives; we use a 
fermion field with Yukawa interaction to an external field as an example which is 
sufficiently simple and sufficiently general. We compare the computed functional 
derivative to analytical estimates. 
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1 Introduction 



The one-loop effective action plays an important role in various contexts of quan- 
tum field theory. It leads to a correction to the energy of classical solutions (see 
e.g. IH) and to corrections to quantum tunneling and thermal transition rates 
[3-5]; it is the basis of Hartree type selfconsistent field equations. If in the for- 
mer case the corrections are important the effective action should be included 
in determining the extremal classical field configuration. In Hartree type com- 
putations as e.g. in the nontopological soliton model of Friedberg and Lee ||, 
or in the Nambu-Jona-Lasinio model [8-11] the classical action may even have 
no solution and the classical field may have no genuine action, it may be just a 
parametrization of condensates. Denoting the "classical" field as (j) the equation 
of motion reads 

5T 5Sr, SS}-^ 



^1 +^ = (1.1) 



50(x)50(x) 50(x 
where Sd may be absent. The one- loop effective action takes usually the form 

, 1, det(-A + m^ + \/(0(^))) 
^'^^ ~2^'' det(-A + m2) ^^'^^ 

where the sign + or — refers to bosonic and fermionic actions respectively. The 
Klein-Gordon type operators will in general be multichannel operators, mass term 
and potential being given by matrices. Two methods for computing the effective 
action numerically have been proposed recently in [12,13] and in [14-19] and have 
been applied in various contexts [20-24]. Both methods are closely related. In 
both of them the usual problems of renormalization and of removing zero and 
unstable modes are easy to handle. While the second method , based on a theorem 



on functional determinants ||25[, is considerably faster the first of them is more 
convenient in the present context. It can be adapted easily for computing the 
functional derivative 5Slg^/6(f){'x). Renormalization remains easy to handle and 
will be discussed. We will not consider zero or unstable modes here. These are 
particular to specific models and here we would like to keep the presentation of 
the method general. For reasons of economy we restrict ourselves furthermore to 
the case of one scalar external field; this restriction can easily be removed. For 
the fluctuating field we choose a Dirac field as the simplest case of a coupled 
channel problem for the fluctuation operator. 



2 The fermion determinant and its functional 
derivative 

We consider a fermion field in 4 Euclidean dimensions which receives its mass 
irtp = gv from its coupling to a Higgs field 0(x) which takes the vacuum expec- 
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tation value v. The Dirac operator has the form 

V = ^- g^{x) = I?W - - v), (2.1) 

where V^^^ = ^—mp is the free Dirac operator in the 'broken symmetry phase'. V 
is the Euclidean Dirac operator, the gamma matrices are hermitean and satisfy 
{lfj.,lA = 25^,^. Since V is not a positive definite operator we introduce the 
operator 

M=j5'D-f5'D = -d^ + ml + V{x), (2.2) 

with the potential 

Vix)=g^i^\x)~v^)+gmx) ■ (2-3) 
The Euclidean efi^ective action induced by the fermion field is then obtained as Q 

SM = -lndetr^U-^lndet^-^^+-- + ^^^ 



= -^Indet (l + (-9^ +m|)~V(x)) . (2.4) 

As a next step we will derive an explicit expression for the functional deriva- 
tive of 5'efr[0] with respect to 4>{x). As a starting point we use the perturbative 
expansion of the effective action with respect to the potential V{x). We define 
the free Green function Qo{x — x') as the solution of the equation 

i-d^ + rriDGoix - x') = {x - x') (2.5) 

where is the Euclidean Laplace operator in 4 dimensions. Then the effective 
action can be expanded as 



oo / -i\n „ n 

Sea[<l)] = ~^ — J n '5o{xi - Xi_i) V{xi) (with xq := Xn) ■ (2.6) 

n=l 1=1 



Obviously only V{x) depends on 0(x), so the functional derivative is obtained in 
a straightforward way as: 



6S,s ^ (-1)\ f,^ 5V(xi; 

— " - j \ ii'^xiyoyxi ^ -^i-i) '^v-^i) j yoyxi ~ 

(2.7) 

If one looks for static solutions, the actual variation will be with respect to fields 
depending only on the three dimensional vector x. We will restrict our consid- 
erations to this case in the following. For instanton computations (see e.g. 
one introduces an extra Euclidean time variable, so that the present restriction 
represents no loss of generality. Obviously the time integrations can be carried 



^We omit the superscript 1 — 1 in the foUowing. 
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out explicitly. In order to do so we introduce the Green function ^o(x — x', i/) 
which is related to Qo{x — x') via a Fourier transform with respect to Euclidean 
time T. It satisfies 

( A - m| - i/')^?o(x - x', u) = (x - x') (2.8) 
and is given explicitly as 

/X /■ d-Vi exp(ipi(x-x')) . . 

- " ' = y (2^ pf+.^ ^'-'^ 

with := m\ + i/^. The time integration leads to delta functions for the differ- 
ences of the Euclidean energies i^j, leaving just one integration with respect to v. 
Explicitly we find 

= 27r5(0)f-i)tr/d^x,^ (2.10) 
(5(/)(x) ^ ' \ 2) J 5(f){:x.) ^ ' 

/ ^ E(-ir^Mnd'^^^o(xi-x,_i,i.)V(x,)Uo(xi-x„,i/). 

_"L n=l \i=2 / 



oo 

X 

-oo 



The factor 27r5(0) is due to a delta function 27r5(z/„ — z/q) whose argument vanishes 
due to the cyclic identification of Vn with vq. It is to be identified, as usual, with 
Euclidean time via 27r(5(0) — > r. Indeed, for time independent external potentials 
the effective action is proportional to r and in fact we are considering Ses/r, the 
zero point energy. As Gq{'x — x', i/) is an even function of u the integral over this 
variable may be replaced as 

oo oo 

J du — >2jdiy. (2.11) 

-oo 

Furthermore the perturbative series in the potential V(x) can be recollected into 
an exact Green function 

oo / „ n \ 

6;(x,x',z/)E(-l)"+i / nd^a;,^o(x.-x,_i,z/)V(x,) k;o(x'-x„z/) (2.12) 

n=l V"^ i=2 J 

with Xi = X. ^(x, x', u) satisfies 

(A-/«2-V(x))6?(x,x',i/) = -5(3)(x-x'). (2.13) 
With these replacements the derivative of the effective action becomes 
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It remains to evaluate the derivative of the potential V(x) with respect to 0. One 
obtains 

= (2gV(xi)+97y„)^<3>(xi-x). 
Inserting this expression into (|2.14|) yields the result 



1 SSes g 



2, 



r 6(f){x.) TT 



oo oo 

/ dz/tr^(x,x, z/) + — / dz/tr7V^(x,x, z/). (2.15) 

J 271 J 



Here the gradient in the last term acts on the Green function at equal arguments, 
i.e. one takes the gradient after taking the limit of equal arguments. The expres- 
sion we have obtained is a formal one, yet. The Green functions are not defined 
at equal arguments, they have to be regularized and a renormalization has to 
be performed. This will be discussed below; we prepare this here, however, by 
introducing a convenient notation. As we have seen ^(x, x', u) can be expanded 
with respect to V(x) via (|2.12| ). The trace with respect to x and with respect to 
internal indices is equivalent to a summation of one loop Feynman graphs. In the 
following we denote by a superscript (k) a contribution of order k in the poten- 
tial V(x) and by the superscript (k) the sum of all such contributions starting 
with order k. With this convention we can define a finite part of the functional 
derivative via 

1 6S^ ^V(x) 7,,_^M.,. „ , 9 ^ 



T 50(x) TT 



y dz/tr^(2)(x,x,z/) + ^ y dz/tr7V^(2)(x,x,z/). (2.16) 



This quantity can be evaluated numerically; the first and second order contribu- 
tions are divergent and will be considered later. 

In some applications the effective action is computed in a quantum state in 
which some of the fluctuations are excited, e. g. in the form of occupied valence 
quark states. Then the "sea quark" Green function has to be modified by adding 
the "valence quark" contributions 



:e™'(x,x,i/) = ^_ V,„(x)^.t(x) (2.17) 



where a runs over the occupied levels. Since the operator 7V — g<f){^) occuring 
in ( p.l5[ ) is just the static Dirac operator and 

(-7V + (70(x)) ^,(x)E„7o^„(x) (2.18) 

this means that we obtain an additional valence contribution 

r50(x) ^ 

as to be expected. 
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3 Computation of the Green function 



The computation of the Green function ^(x, x', u) has been considered previously 
in Ref. [0, |r^. In order to proceed we have to introduce a further restriction 
on the classical configuration 0(x). We have to assume that it is spherically 
symmetric. For most cases of interest this is the minimal or extremal 

solutions have usually the highest possible symmetry. Then we may expand the 
exact Green function into partial wave Green functions. These can be obtained 
by solving the radial Schrodinger equation numerically. 

3.1 Partial wave reduction 

For the partial wave reduction we have to use an appropriate spherical spinor 



basis. It can be found in |]26|. One introduces the 2-spinors 

forj = / + l/2, (3.1) 



_ ( V^ + l/2 + mr-/^(x) 



= -iK^i 'rv^^^r^^.m-,:;'. I for J = /- 1/2.(3.2) 



I / ^/ + 1/2 -mY;^'^^\x) 

which are eigenspinors of J^,Jz, and S^. Further they satisfy 

ax(^(±) = (3.3) 

Using these spinors one can construct 4-spinors and reduce the Dirac equation 
into two coupled radial equations for the partial waves. In order to obtain a 
similar Ansatz for the Green function we introduce the matrices Qj^\x,x') via 

nf{x,x'):= t vf2{x)^f2\x'\ (3.4) 

m=-j 

in terms of which the Green function can be expanded as 

g{±)l^_^ x' i^) - ^ f ^iS'^'^' ^f\^^ -9ti2{r, r', u) nf\x, x')iax' 

3 \ 9jfi {r, r', u) iaxnf^ (£, x') (r, r', u) il^p {x,x') 

(3.5) 

The partial wave Green functions Qj^mi^^ ^'^ ^) satisfy 

((5„fc(-D„ + «:2) + Vnk{r)) gi^ir, r', v) = S^ J^""/^ . (3.6) 



Here Dn is the radial operator 



= ^ + 2 A _ + 1) ,3 
" dr"^ r dr 



5 



with /i = J - 1/2 and /a = j + 1/2 for and /i = j + 1/2 and /s = j - 1/2 for 
(7j \ The potential has been reduced to a (2 x 2) matrix 

One sees easily that the (+) and (— ) systems at equal j are equivalent upon 
identifying 

9ll 922 ; 922 9ll ; 9l2 ^ 921 ; 921 ^ 9l2 ■ ['^■^) 

3.2 Partial wave Green functions 

The partial wave Green function can be expressed in standard way by regular 
and singular solutions of the homogeneous equation 

E [^nki-Dn + K^) + Vnkir)] Mr) = 0. (3.10) 

k 

This system has a fundamental set of two independent two-component solutions 
regular as r ^ cxd which will be denoted as Z^"*" (a = 1, 2) and two such solutions 
regular at r = 0, denoted as . In terms of these solutions the partial wave 
Green function can be expressed as |jl3| 



(3.11) 

where ojap is related to the Wronskians of the two sets of solutions via 



E {fn^i:fn" - ft^^fr) = r'WirHr), f-{r)) . (3.12) 



A short derivation of this formula is given in the Appendix. It is convenient to 
choose the two sets of solutions in the form 

f:^ = {s: + h:^{r))bi{^r). (3.13) 

Here bf{z) = ki{z) and bY{z) = ii{z) are the modified spherical Bessel functions 
regular for z oo and as z — > 0, respectively. The functions /i"^ are regular both 
at r = and as r oo. They can be chosen to vanish as r ^ oo as a boundary 
and normalization condition. Then the Wronskian simplifies and we have 

Ucfj = 5a(3/l^ . (3-14) 



''We suppress the angular momentum index j, the superscript (±) and the argument v in 
the foUowing. 
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We will need the Green function only at r = r'. The limits r' \ r and r' r 
are equal by construction, they involve however different sets of functions. The 
numerical differences are found to be unimportant. We have chosen to take the 
symmetric limit 

gmn{T.r) =Y.^ {/rW/r W + /r (r)/r(r)} . (3.15) 

3.3 Removing the leading orders in the external potential 

As mentioned in section 2 we will have to remove the leading orders in the external 

potential V(x) in order to obtain the finite part bSf^ j b(\){y:)\ this implies obtaining 
the Green function in order (2). We note that our ansatz ( |3.13| ) for the solutions 
/"^ represents a separation into the zeroth order part proportional to 5° and the 
order (1) contribution proportional to the functions . These functions satisfy 
the coupled differential equations 



(I. + 2 + j I) KHr) - E (C . ^ . (3.16) 
which can be written in symbolic form as 

m^ = y(i + /i±). (3.17) 

If one expands the functions h with respect to powers of V as 

:= h^~) = h^W + h±i2) + . . . ^ (3.18) 

one obtains the following sequence of equations 

p/i±(i) = V {l + h^'^) , (3.19) 

= V, (3.20) 

= y/,±a) (3.21) 

which can be solved numerically. In terms of these functions the partial wave 
Green function at r = r' can be expressed as 

9nm{r,r) (3.22) 



The third and fourth term represent the contribution of order (2) 



(3.24) 



3.4 The derivative of the effective action 



Having at our disposal a numerical method to compute the Green function it is 
straightforward, using ( p.l6|) , to derive an expression which allows the numerical 
computation of the functional derivative of the effective action. We first have to 
take some traces of the Green function and its derivatives. Using the identities 



J + 1/2 / 1 



47r 



1 



(3.25) 



and 



I 



m=—l 



m\2 



21 + 1 
= 



j 



ll + -~mJl + - + mYi 



m+l/2 ^m-1/2* 
1] 



one finds 



-| oo 

trG'(2)(x,x,z/) = - E (2J + 1) 

4^ .=1/2 



1 oo 
i=l/2 



(3.26) 
.(3.27) 



Supplying a factor of 2 in order to take into account both the / = j ' + 1/2 and 
the / = j — 1/2 modes we obtain the final result 



r 50(x) 



^V(^) 
27r2 



dz/ E (2J + 1) 11 (^,0+^22(^,0 

i=i/2 



27r2 



dz/ E (2J + 1) ^7;S(r,r) + -(?g^2(r,r) 
i=i/2 L r 



(3.28) 



This expression is finite and can be computed numerically. However, we have 
to include the renormalized divergent parts as well. This will be done in the 
next section. Here we will discuss some details of the numerical evaluation and 
we present an analytic approximation, based on the derivative expansion. This 
analytic approximation will be compared to the numerical results. 
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3.5 Some details of the numerical evaluation 



The numerical evaluation of the effective action has been described in extenso 
in previous publications ||T2| , [13| . If one computes the effective action the first 
operation is to integrate the trace of the Green function over the radial variable 
r. Thereby one obtains the partial wave contributions at fixed v which are then 
summed using an asymptotic extrapolation. Finally one has to do the integration 
over V. Here the information we want to draw from the Green function is a local 
one, we want to determine 5S'cfr/(50(x) as a function of r. This requires a different 
order of steps. Let us write (|3.28| ) in the form 



where now 



1 bS. 



(3) 
eff 



r 50(x) 
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27r2 



p OO 

y dz/ ^ (2j + l)t,(z/,r) 



(3.29) 



i=i/2 



+ 



'(2) / \ 2 (2) , ^ 



(3.30) 



The fact that the method for computing the effective action developed in |T^, [13 



is very fast is due to the fact that after only two Runge-Kutte runs for computing 
the functions /i^^(r) one knows the partial wave Green function at all mesh points 
r„. This advantage can be maintained in the following way: the Green function 
is computed, at fixed angular momentum j for all r„ in one sweep and at each of 
the (roughly 2000) mesh points r„ the appropriate traces tj{i',r) are computed 
and stored as functions of u and j. 

The computation was extended to j — 1/2 < 25 and to z/ < 2. The partial wave 
sums 0"„(z/) = Z](2j + l)tj(z^, Tn) were computed using an extrapolation of the form 
a/j^ + b/j^ to include the terms with j — 1/2 > 25 and the integration over u 
was extended beyond z/ = 2 by using an extrapolation of the form a/z/^ + 6/z/^. 
The form of the extrapolations can be derived rigorously, the convergence was 
checked for several values of r. 

The total time for one evaluation of SS^^^ /5ip{x) at all values of |x| takes 
about 5 minutes on a workstation IBM Rise RS6000/340. 



3.6 Analytic approximation 

In |T5[ we presented an analytic approximation to the finite part S^g of the action, 
using an expansion with respect to gradients of the external field 0. It reads 



9^ 



IGtt^ 

-.2 



In- 



2\2 



9 



167r2 



V 



3v^ 



(3.31) 
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This approximation was found to be excellent at sufficiently high fermion masses. 
It is interesting to compare exact results and this analytic approximation also on 
the level of the functional derivative. Taking the fields to be time independent 
and dependent only on r = |x| we obtain 

-v')^ . (3.32) 

This analytic approximation is compared to the numerical results in Figs. 1 to 
4. For the numerical evaluation we relate the physical quantities to dimensionles 
ones hy X = x/v and (j) = vcf). Here and in the following we use for the profile 
function 

/(r) = (l + exp(-y;^2^+y:R2^))-i (3.33) 

with i? = 2 as in [|15|. It has the typical shape of a medium-sized bubble. The 
approximation is found to reproduce the numerical results quite well; as to be 
expected the approximation improves with increasing fermion mass. 

To the finite part of order (3) we have to add the finite parts of the first 
and second order contributions. They are known analytically and their numerical 
evaluation, involving only a numerical Fourier transform and simple integrations, 
can be performed with high precision. For this reason we have not included them 
here into the comparison between numerical and approximate analytical results. 



1 5S. 



(3) 
eflf 



9 



T 50(x) 



47r2 
- 2 



'^In^- 

V 

r 



167r2 



In — 



2 

3^ 



4 Renormalization of the effective action 

Having discussed the computation of the finite contributions of third and higher 
order in the external potential V(x) we have to come back now to the compu- 
tation of the first and second order contributions to the effective action. These 
contributions are divergent. As a first step they have to be regularized, e. g. by 
dimensional or Pauli-Villars regularization, introducing thereby an ultraviolett 
cutoff A. If there is no genuine action for the field 0, as e.g. in Nambu-Jona- 
Lasinio type models where (f) parametrizes fermion condensates, the theory will 
depend on this cutoff - the latter will be a physically significant quantity. If the 
(p field has a genuine action of the 0^ type then the divergences can be absorbed 
into counterterms and the cutoff dependence disappears via renormalization of 
the physical masses and couplings. We will consider the latter case here; to be 
specific we choose a Higgs type potential for the field (j). 
The bare Lagrangian of a Higgs- Yukawa theory reads 

L = U^ct>d^'ct> - U{ct>) gct>) ^ . (4.1) 
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with 

f/(0) = ^ (02 - . (4.2) 

If the Higgs field is assumed to take its vacuum expectation value v at large r it 
is convenient to introduce the field shift 

(j){x) = V + ip{x) (4.3) 

and to express the Lagrangian in terms of ip as 

C^d^ipd'^ip - ^ml^p^ - Xvip^ -^V^ + i^ {i^ - {mp + g^)) ^. (4.4) 

Here we have introduced the boson mass = \/2Xv and the fermion mass 
mp = gv. 

The counterterm Lagrangian will have to include a fourth order polynomial in 
the field and also a wave function renormalization counterterm. We do not 
have to include here corresponding counterterms for the fermion part of the ac- 
tion. Such terms would be made necessary only if boson loops were included. The 
counterterm Lagrangian reads then 

C,.,XV) = -A^ - \b^' - -^C^' - ^^D^' + Uz d,^d^^ (4.5) 

For the Yukawa - Higgs Lagrangian we note that the bare Lagrangian depends 
only on the sqare of the unshifted field (f? and its derivative. This holds also true 
for the divergent terms of the one-loop effective action generated by the fermion 
field, due to the structure ( [4.2| ) of the external potential V(x). Therefore, only two 
of the four constants A,B,C and D are independent. This implies the relations 

A-Bv + ^Dv"^ = 

C = vD. (4.6) 

The choice of the independent quantities depends on the renormalization condi- 
tions which will be given below. 

The first and second order contributions to the one-loop effective action are given 
- first as formal definitions and then in dimensionally regularized form - by 

(47r)^ [ e rrip J J 
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1 



647r2 
+ 2- 



1 tJ+p(1) 



i + V'(l)+lni^ 

e nip 



(4.7) 



p{(l) j ' 



Here p{q) is defined as p(g) := q/ \J cf + 4:m'jp. We choose the renormahzation con- 
ditions in analogy to Coleman und Weinberg . In order to do so we introduce 
at first the one-loop effective potential. It is given essentially by the one-loop 
action induced by constant fields ip{x) = (f. Then the external potential takes 
the form 











1 
1 



with 



(4.8) 
(4.9) 
(4.10) 



and the effective action can be written as 

where ^2*^^^ is the volume of space-time. The renormalized effective potential in- 
cluding the bare potential and the counterterms is given by 

Kfr,ren : f/(0 = + ^) + V^^' + {^^ + ^B^' + Icf' + ^D^'} ■ (4.11) 

Here V^g~^ includes the one loop contributions to all orders, i. e. the regularized 
divergent parts of first and second order and the (3) contribution of all higher 
orders. The first and second order terms can be found easily from (^4.7|) ; they read 



1-1(1) 



eff,reg 



m. 



and 



1-1(2) 
eff.reg 



(47r)2 
(47r)2 \ 



- + ^(l) + l + ln 



47r/i2 



m. 



m. 



(4.12) 



(4.13) 



The (3) contribution has been given in |]T5|, it reads 



1-1(3) 



eff 



E 



\n+l 



n 



In 1 + 



2mj, 
(4^ 



m. 



' m y 

^p"^ + nip ) 

1 : 

2 ml 



1V%) 
2 m% 



(4.14) 



2ml 



4 m'p 
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We require the effective potential to retain the following properties of the tree- 
level potential U{(f) = v + (f) : 



dV;fr,ren(<^) 



dip 



= (physical mass), (4-15) 

<fi=0 

= ((^) = 0). (4.16) 

<fi=0 



The first of these conditions fixes the constant B of the counterterm action, the 
second one the constant A. The remaining constants are fixed by the constraints 
([4.6|) . One finds after some algebra 

27r^ \e rrip J 

D = - + V'(1) +ln- ^ * 



With these counterterms the total effective potential becomes 

W„(.) ^ ,4.17) 

-[gip + mF) \n 



(47r)^ 

and is independent of e and /x. In order to determine the wave function renor- 
malization we have to consider again. Inserting the explicit expression for 



tr|(V)(g))|2 into (O) we find 



where V^(y') has been defined in ( |4.9| ). In the integrand the contribution 



\Vmq)?{- + m+^^^] (4.19) 



mi 
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has been absorbed already into the renormahzed effective potential, (see ( |4.13| ) 
and ( [4.17|) ). The contribution 



OTT^ J {2tt)* \e nip J 



167 



(4.20) 



, +V;(l) + ln^ d^x{d,<P{x)y 



will be absorbed into the wave function renormalization 6Z (see below). This 
leaves a finite part 

I p{q) 1 - p{q) I 

The expression in curly brackets can be expanded as 



Therefore the term g'^q'^\(f){q)\'^ in the first parenthesis will be multiplied at least 
with and yield terms with at least four derivatives of the fields. On the other 

hand the term |(t^(0)(5')P will yield another second derivative term when multi- 
plied with q'^/Qnip. We denote this contribution as S^^^^^^-^^^-^^. It can be rewritten 



as 



'S'ifT,analytic — TcU). (4.23) 



levr^ 



Since we renormalize at = w, i. e. <y9 = we find here another contribution to 
the wave function renormalization : 



9 



2 



247r2 



d''x{d^^Y , (4.25) 



/ON 

to be absorbed into 5Z. Other second derivative terms appear in S),^ given in 
( p.31| ), their coefficient vanishes at = f . We find therefore, collecting the con- 
tributions (|]20D and { WM ) 

^Z = -^fi+V.(l)+ta^^)+T#i. (4.26) 

(47r)2 \e ml J Uir^ 
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The effective action 

T = S,i + Sis' (4.27) 

becomes therefore 

While taking the functional derivative of the first term is trivial and the numerical 
computation of the term S^g has been discribed in section 2, we have to discuss, 

(2) 

as a final step, the computation of the functional derivative of Sgg g^j^g. Restricting 
ourselves again to time-independent configurations we find, using the definition 

2/0(x)(02 - vYiq) + q'i<P- ^)^(q)) ^ in ^^44 ■ (4-29) 



p{q) 1 - p{q) 

Here (..)'" denotes the Fourier transform of the expression in parentheses. Taking 
into account of the fact that the fields depend only on r = |x| and therefore their 
Fourier transforms only on g = |q|, this expression can be further simplified. 
Its computation involves simple one- dimensional Fourier-Bessel transforms and 
integrations. 

We compare the numerical results for the renormalized second order contribu- 
tion (l/r)(5S'gfinite/'50(x) -^2A0/127r2 with those for {1/t)5S'^ /6(f) in Figure 5. 
The renormalized second order contribution is found to be considerably smaller 
than the contributions of order (3). 



5 Conclusions 

We have presented here a method of computing the functional derivative of an 
exact effective action F for spherically symmetric background field configurations. 
The vanishing of this functional derivative determines field configurations that 
minimize the action or, for time independent configurations, the energy of a 
quantum system in one-loop order. We have considered the exact computation 
of the finite part as well as the implementation of renormalization, whenever it 
is necessary. We have considered here a fermionic loop contribution, the method 
is of course more general. 

For the special case considered here we have also presented an analytic esti- 
mate 



1 5S. 



eflf 



ml 



T 



02 (in 02 _ 


+ 1 
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^" + ^0') (in 02 -^(0^-1))] (5.30) 

for the finite part of order (3) and compared it to the exact numerical results. 
The approximation is found to be good at large rriF, as to expected. This local 
expression, as well as some local terms occuring in SS^q^^^^^^/ 50(x) (see ( [4.29|) ) can 
be introduced explicitly into the classical equation of motion ( |1 . 1| ) thus obtaining 
an improved differential equation for the classical solution. Its solution should 
represent a good starting point for an iterative procedure which takes into account 
the remaining parts that have to be computed numerically. 
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Appendix 



Equations and (|3.12|) are relations one expects to find in some textbook, 



therefore in [jIB|| they were given without proof. Since we have not found such a 
textbook yet, we give here a heuristic derivation. The Green function we want to 
determine satisfies 



5nki-^, --^ + M^!if^ + k') + v^,{r) ] gUr. r') = -^5{r - r') . (A.l) 
dr^ r ar r'^ / 



From the the fact that it is a solution of the homogeneous differential equation 
for r 7^ r', both as a function of r and r', and imposing regularity at r = and 
oo we conclude that it must have the form 

gUr. r') = E {e(r - r')f^+{r)c^^ft{r') + Q{r' - r)/^ (r)^.^/^ (r')} 

(A.2) 

in terms of the functions fn^{r) defined in section 3.2 . The Green function has 
to be continuous at r = r'. Its first derivative has to have a discontinuity so that 
the second derivative produces the appropriate 6 function. These two conditions 
lead to the relations 



E(ca/3/r(0/n^-(r)-rf«/3/r(r)/r(r)) = (A.3) 

r'EMrw/r(r)-rf./,/rw/r(r)) = s^n. (a.4) 



a,l3 



We contract the first of these equations with r'^f^^ir) and the second one with 
f2i^{r) , then we subtract the resulting equations from each other. We obtain 

E {c^.r'Wip-', r^)ft - d^pr'WiP^, r-)fT{r)) = -/„^+(r) (A.5) 

where we have defined the Wronskians as 

W{r, /) = E {fZMir) - C^'-^giir)) . (A.6) 

m 

If the potential Vnk{f) is symmetric the differential equation satisfied by the func- 
tions f^^ir) implies that r'^W{f", g^) is independent of r. Since the functions 
fn^{r) vanish as r — > cxo the Wronskians W{f'^^ , f""^) vanish. Using the fact 
that the solutions f'^'^{r) have been chosen to form a linearly independent set, 
we conclude that 

E d.pr'w{r+, r-) = E ^-z^^"" = (a.t) 



17 



or 

dap = ^al ■ (A.8) 

In order to obtain the coefficients Cap we contract (A. 3) with and (A. 4) 
with f^{r). Using the fact that vanishes because it does so as 
r — i> 0, we find 

Ca/3 = ^^pa ■ (A.9) 
This closes the demonstration of (|3.11|) and ( p.l2|) . 
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Figure Captions 



Fig. 1 The functional derivative of the finite part of the effective action (l/r)(5S'gg /50(x) 
in units v^, calculated for a spherically symmetric time independent profile 0(x) = 
/(r) (see (p.33|) ) with R = 2. Here we choose the fermion mass in units of v to 
be rhp = g = 1.5. The solid line is the exact numerical result and the dashed line 
is the analytic estimate given in Eq. ( p.32| ). 

Fig. 2 The same as in Figure 1 for = 1.7. 

Fig. 3 The same as in Figure 1 for (7 = 2.0. 

Fig. 4 The same as in Figure 1 for g = 2.2. 

Fig. 5 The various contributions to the full one-loop correction of the functional 

derivative of the effective action in units of for q = 2.0: The dotted line 

(2) 

represents the renormalized second order contribution (1 /T)5S^ff ^^^^^^f 5(f){x) — 

(7^A0/127r^; the dashed line is the exact numerical result of {l/T)6S^g /6(f) which 
is also presented in Figure 3. The solid line is the sum of both contributions. 
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